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The process of top quark pair production at the ILC with subsequent decays of one of the top quarks to a charged 
Higgs boson and b-quark is considered. The charged Higgs decays to tau leptons whose polarization is the opposite to 
those coming from W bosons. This difference is reflected in the energy distributions of the tau decay products in the 
top quark rest frame, which can be reconstructed at the ILC using the recoil mass technique. We present an analysis 
including spin correlations, backgrounds, ISR/FSR and beamstrahlung, and show that a fit of the shape of the pion 
energy spectrum yields the charged Higgs boson mass with an accuracy of about 1 GeV. 
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1. CHARGED HIGGS AND TAU POLARIZATION 

The mechanism of spontaneous electroweak symmetry breaking in the MSSM leads to a Higgs sector with five 
physical states: two CP-even Higgs bosons, h and H, one CP-odd one, A, and two charged Higgs bosons (If CP 
is violated then the three neutral Higgs bosons will not have definite CP-parity.). The discovery of the charged Higgs 
bosons would be very important, as it would show directly that the Higgs sector has more complicated structure 
than the one in the Standard Model, thus providing clear evidence for physics beyond the SM. 

The charged Higgs boson couples strongly to the fermions of the third generation. If the charged Higgs is heavy 
{Mfj± > Mt), it can easily be detected in decays to top and bottom quarks {H^ tb) by examining the tb invariant 
mass spectrum. In this case, Mh± may be measured to, at best, about 5 GeV at the LHC [1], and about 1 GeV at 
an e'^e~ linear collider running at an energy of 800 GeV and collecting a luminosity of 500 fb~^ [2]. In a number 
of other scenarios, however, the charged Higgs bosons will be rather light {M[j± < Mt). A precise measurement of 
Mh± in that range is a challenging task for any collider. In this case the charged Higgs decays dominantly to a tau 
lepton and neutrino {H^ t^i^), and it is impossible to reconstruct directly the invariant mass of the di-tau final 
state. However, due to the polarization of the t"^ leptons, the energy of the t"^ decay products depends strongly 
on Mh± , a feature that can be exploited to extract Mh± indirectly. The main background to a signal in the 

decay mode comes from the VF-boson decays t^v. However, thanks to the different structure of the 

and electroweak interactions to r^-leptons, the tau polarization is very different which will allow a separation 
of iJ± T^v and r± on a statistical basis. More specifically, the decay products (r+ — > 7r+i^, or 

T+ — !■ p^V^ etc.) have strikingly different topologies according as to whether they originate from a parent or 
. The importance of tau polarization in searches for charged Higgs bosons has already been stressed in previous 
studies [3]. 

The key point is that leptons arising from t^v decays are almost purely right-handed, in contrast to 

the left-handed leptons which arise from W~ decays. This contrast follows from the helicity fiip nature of the 
Yukawa couplings of the Higgs fields, and the helicity conserving nature of the gauge interactions. The most dramatic 
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Figure 1: tt meson energy spectrum in the resonance rest frame. 
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Figure 2: The energy spectrum of the vr* meson (left) and meson (right). The dotted hne corresponds to the background, 
and the sohd one to signal. 

difference is seen in the energy distribution for the single pion channel {H^ /W^ —> t^v — > tt^vv) in the rest frame 
of the parent boson (W^ or H^) - see Fig. 1. 

In practice, it is nearly impossible to reconstruct the rest frame of the or bosons because the momentum 
of the neutrinos cannot be measured. Instead, one can use the top quark pair production channel with top quark 
decays {t — > H^b) and reconstruct the top quark rest frame. This is possible at a linear collider when the second 
produced top quark decays hadronicly, in which case one uses the recoil mass technique analogous to those in Higgs 
boson studies [2, 4]. 

The distributions for the energy of the single pion in the top rest frame has the following form [5] for and 
cases, respectively. 

r dy^ X„rax - X„»n \(1 - p^)log^ + 2P,(1 - ^), X„,„, < 

where y„ = a;™„, = ^jf^, x^ax = -ET™ = Ssf^' = '^'^^ t^*^ ^ boson, P.^ -1, and for 

the charged Higgs boson, P^- = 1. ^ 

The energy distribution for a pion from decay has a maximum at the point E{tt^) = M//±^/(2Mt), as shown 
in Fig. 2 (left). This dependence allows one to extract Af^± from the shape of the spectrum. One can also take into 



^For simplicity, in Eq. 1 we liave neglected tlie b-quark and tau-lepton masses, while these masses liave been included in all numerical 
simulations. 
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account decay channel (r+ — > p'^y), which has twice the branching ratio as the pion channel. The shape of the 
energy distribution is more complicated, however, and is less sensitive to Mjj±, as shown in Fig. 2 (right). More 
detailed studies are needed for this channel. 



2. EFFECTIVE tbH^ INTERACTION IN THE MSSM 

For our simulations described in the next section, we use the effective Lagrangian approach for the MSSM presented 
in Ref . [6] . The effective Lagrangian of charged Higgs interaction with fermions of the third generation has the form: 

where fhh is the running bottom mass in the MS scheme. The width of the top quark decay to a charged Higgs and 
b-quark takes the form: 

T^ssMit^bH-)^ Tl + Am.)^ (3) 

where 

Am, = ^ tan/3 ~ ^^^^V^ tan/J + (4) 
hb on max{mi ,mt ,M~) ° 

0102 " 

The Amfo correction proceeds from the one-loop vertex corrections, which modify the relation of the bottom Yukawa 
coupling to the bottom quark mass. Similar corrections are induced in the neutral Higgs sector [7]. The QCD 
improved top-quark decay width is given by [6], 



1 + 



7 - — - 21og(l - + 2(1 - 
+ (^ + ^log(l-g//+)J (l-te+f 



(6) 



The above expression, Eq. (3), takes into account the dominant supersymmetric corrections to all orders in per- 
turbation theory. It also includes the rcsummation of the large mt/mf, and m^/m;, logarithms, which amounts to 
replacing rrib by the running bottom mass at the proper scale [6] , as well as all finite one-loop QCD corrections. 

In the numerical similations, for any given set of supersymmetric particle mass parameters, we have computed 
the Amfo corrected charged Higgs Yukawa couplings with fermions of the third generation, by means of the program 
CPSuperH [8]. 



3. SIMULATIONS 

The numerical simulations have been performed assuming a center of mass energy ^/s = 500 GeV and a total 
integrated luminosity C = 500 fb~^. We have performed detailed computations and Monte Carlo simulations for 
three different sets of MSSM parameters, all leading to light charged Higgs bosons. The first two scenarios are based 
on the mass parameters Mq = Mu = Md = I TcV, Mj = M2 = I TcV, At = 500, /i = ±500 GcV, and tan (3 = 50, 
which give Mh± = 130 GeV. The branching ratio BR{t H'^h) will be enhanced or suppressed depending on the 
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sign of /i. The third parameter set falls in the so-called "intense-coupling regime" where the neutral Higgs boson 

masses are close to each other, and they couple strongly to fermions of the third generation [9]. In this case, we have 
jjL = 1000 GeV, tan /3 = 30 and Mh± — 146 GeV. The branching ratios for these three examples are 



(i) M//± = 130 GeV, /X < and tan /3 = 50 
(a) Mh± = 130 GeV, > and tan /? = 50 
(in) Mh± = 146 GeV, a* > and tan /3 = 30 



BR{t H+b) = 0.24 
BR{t H+b) = 0.091 
BR{t H+b) = 0.063 



These points are not excluded by Tevatron searches which lead to the bound BR{t H+b) < 0.42 at 95% C.L. when 
Mh± < 150 GeV [10]. The couplings have been implemented in CompHEP [11], which has been used to compute 
the cross sections for signal and background processes, including decays of top to W+ and H+ which subsequently 
decay to polarized t+ leptons. CompHEP was also used to generate events, and effects from initial state radiation 
and Beamstrahlung were included. Polarized t+ decays have been simulated using TAUOLA [12] interfaced to 
CompHEP. Hadronization and energy smearing in the final state are accounted for by means of PYTHIA [13] using 
the CompHEP-PYTHIA [14] interface based on Les Houches Accord [15]. Effects from final state radiation have 
been implemented using the PHOTOS library [16]. 

A brief description of our fitting procedure is given here. The fitting function is the sum of two functions, one 
for signal and one for background. For simplicity and without loss in accuracy we chose the following form for the 
charged Higgs signal piece, which is motivated by the theoretical spectrum: 

H{x) = pariko{\x — ki\ + - k^), x ^ E^^ (7) 

where: fco = Mf /[parjiparj - Mf)], ki = parl/{2Mt), k^ = 2parl/M'^ - 1, fcs = par|/(2Mf), with par2 = Mh±. 
There are two free parameters: the parameter pari gives the overall normalization of the signal events, and the 
parameter par2 is the charged Higgs mass itself. The function H{x) has zeros at the points E„ = and E„ = Mt/2, 
and a maximum at the point E^^ = Mfj±/{2Mt) as follows from the theoretical spectrum. For the background 
function we use a forth-order polynomial of the following form: 

B{x) = parsfix), x = E„ (8) 

where f{x)—co+cix + c^x^ + c^x^ + 040;^, and the fitting parameter par 3, gives the normalization of the number of 
the background events. The coefficients cq, . . . , C4 are fitted from the shape the background pion energy distribution. 

We use the method of maximum likelihood to fit a spectrum created from simulated signal and background events, 
and obtained the following results for the three cases described above ^ 

(i) Mh± = 129.7 ±0.5 GeV, 
(n) Mh± = 129.4 ±0.9 GeV, 
(in) Mh± = 145.5 ± 0.9 GeV . 

Figure 3 shows results of the simulations and fits. 



4. CONCLUDING REMARKS 

The tau polarization is a powerful discriminant separating H+ t+v from W+ t+v in top decays in scenarios 
with rather light charged Higgs bosons, and with large values of tan/3, which enhance the tbH+ coupling. In this 
study, we have considered three cases with Mh± w 130 GeV and tan/3 = 30-50. They illustrate how the tbH+ 
coupling can be enhanced or suppressed to the extent that BR{t H+b) varies from 6% to 24%. A fit to the 



^Wc arc taking into account the expected accuracy, of about 100 MeV, for the top quark mass determination from the top-quark pair 
threshold scan at the ILC (see Ref. [2]), which adds only a small contribution to the charged Higgs mass uncertainties. 
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energy spectrum of the pion in the r — > it v channel aUows one to infer M^± with an uncertainty at the level of 
0.5-1 GeV. This study is a theoretical level analysis, and no systematics or detector effects have been included. On 
the other hand, only one decay mode (r* — > ir^v) has been used, and clearly if one would use other decay modes 
(such as the two-dimensional decay distributions of the p^v and afv decays), the sensitivity to M^± would 

improve. Further work in this direction remains to be done. 
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